Electronic Transmission Through Metallic Nanowires: 
Generalized Scattering Matrix Approach 
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An easy to implement and powerful method for the solution of 3D scattering problems that can 
be well described by Helmholtz equation is presented. The matrix algebra used provides excellent 
stability versus the number of junctions as well as great computational speed. The matrix truncation 
method yields an easy single-parameter convergence procedure. Subsequently, some aspects of the 
electronic transport through metal nanowires are studied by the use of Landauer's scattering ap- 
proach to the conductance. We predict the existence of current vortex-rings patterns due to sharp 
enough narrow-wide connections in atomic size point contacts. Longitudinal resonances between 
scattering centers provide a simple physical picture for the understanding of negative differential 
resistance in ideal monoatomic contacts. Relatively long nanowires with high geometrical perfection 
-like those recently observed by Transmission Electron Microscopy- are modelled exhibiting resonant 
tunnelling and total reflection at given incident energy intervals. 
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Future daily-life technology will be greatly based on 
miniaturization of both electronic components and op- 
tical devices. State-of-the-art research has already pro- 
totyped tools using simultaneously both electrons and 
photons confinement. For instance, in artificially created 
electromagnetic micro-cavities interaction light-matter 
can best be described by new concepts like those of 
"cavity polaritons"EJ, devices in which the light emis- 
sion process can be extra-ordinarily fast. Also, quan- 
tum well lasers utilizing two-dimensional electron gas 
and optical modulators with quantum confinement ef- 
fect are now practically used. The ultimate integration 
of any opto-electronic device has to overcome the theo- 
retical and practical problem of thc^waveguiding at tiny 
sizes for both electrons and photonso. Consequently, it is 
not surprising the huge research effort devoted, in the 
last decade, to these waveguides themselves: metallic 
wires at the nanometer scale and optical waveguides at 
the micrometer scale. A common characteristic lays on 
the basis of these waveguide systems: elastic scattering 
processes are dominant for transport. Metallic atomic- 
size wires and micro-optical guides present lateral dimen- 
sions of the order of the wavelength associated to the 
particles passing through. Helmholtz equation governs 
elastic flux of waves through a general dispersive system 
(note here that the independent free-electron stationary 
Schrodinger equation is equivalent). Therefore, the abil- 
ity to calculate solutions for this equation is of great im- 
portance. 



As mentioned above, the experimental and theoretical 
study of metallic nanowires, namely their electrical and 
structural properties, has been focus of intense research 
in the last decade (see e.g articles in,Ref.0 as well as the 
various experimental works of Reffl, and the theoreti- 
cal works intl). However, only very recently, the com- 
bination of Scanning Tunneling Microscope (STM) with 
Transmission Electron Microscopy (TEM) has .achieved 
atomic-resolution images of metallic nanowiresocj. In 
particular, ultrathin gold wires have been observed ex- 
hibiting a striking geometrical regularity in some cases. 
Here we intend to look at some possible effects due to 
the ideal geometry of these nanowires. Due to the sim- 
ple -therefore general- phenomenology studied, the same 
physics is expected to be found by more sophisticated 
approaches. 

The paper is organized as follows: first we present 
a Generalized Scattering Matrix method for the solu- 
tion of Helmholtz (or Schrodinger) equation in section 
|j} Though the formalism developed can be applied to 
general scalar wave propagation problems and even scat- 
tering of light in some cases, we use it later to address the 
problem of conduction through atomic-size contacts: In 
section [n] we focus on the study of electronic conduction 
through metallic nanowires with several models. Finally, 
in section [ni| we summarize our results. 
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I. SCATTERING MATRIX METHOD 

As a start, let us write the 3D time-independent 
Schrodinger equation in Xp and Ep (Fermi wavelength 
and Fermi energy) units: 

V 2 $(x, y, z) + 4tt 2 [E - V(x, y, z)] $(x, y,z) = (1) 

A general waveguide geometry can be formed by dis- 
cretization in constant section slices (where the confining 
potential V(x, y, z) only depends on x, y). For a quan- 
tum waveguide of constant section, the solution of the 
time-independent Schrodinger equation is separable as: 
<&(x, y, z) = Z(z)ip(x,y), being the separated equations: 

Z" + k 2 Z = (2) 



V 2 ip + [{4ir 2 E- k 2 ) -4ir 2 V{x, y)]ijj = (3) 

where the longitudinal wavevector k is the arbitrary sep- 
aration constant. The general solution for the propagat- 
ing part of the wave function is: Z(z) — Ie tkz + Re~ lkz , 
where / and R are arbitrary constants to be determined. 
In XY we have a generic quantum well determined by the 
confinement potential V{x,y). We will assume that this 
well has some "hard wall" boundary conditions at a given 
limit. The discrete set of eigenvalues {e n }n=o,i,2... (being 
e n = 4:Tr 2 E~k 2 ) can be found analytically or numerically. 
The known solutions of equation [| form a complete ba- 
sis of transversal eigenf unctions A = {<pn{x,y)} n =o.i,2... 
inside the hard- wall limit, so one can write the complete 
solution for the uniform section waveguide as an expan- 
sion in this basis: 



$(x, y,z) = Y J 4>n{x, y) {I n e lk ~ z + R n e- lk ^} (4) 

n 

where the quantity in between brackets can be seen as 
the correspondent complex coefficients of the expansion. 
This leads to a vectorial notation in which the coordi- 
nates of $ in A are given by {I(z) + R(z)}, where the 
j-th elements of vectors / and R are Ije lkjZ and RjC~ lkjZ . 
The basis A is completely characterized by the confining 
potential V(x, y). 



*l-=K A {T A {z)~R A {z)} 

d<S> E 



(6) 



<^=K B {I B (z)-R B (z)} 

where diagonal matrix K A is given by (k A ) 2 = —e A + 
4tt 2 E (analogous definition for K B ). 

Continuity condition for the wave function and its nor- 
mal derivative at the junction leads to these following two 
matrix equations: 



C{B -» A){I B + R B } 



(7) 



C 1 [B -> A)K A {I A - R A } = K {I - R B } (8) 

Here, matrix C characterizes the intermode scattering 
or mode coupling due to the step discontinuity and is 
defined by overlap integrals. The change of basis ma- 
trix to change from basis B to basis A, C(B — > A), can 
only be defined if the eigenvectors of the basis B can be 
expanded in terms of basis A; this implies that the def- 
inition domain for A should include the one for B i.e., 
in geometrical terms, the hard-wall closed line of B is 
completely inside the one for A (see Fig. |l|). 

The element (J, m) of the change of basis matrix 
C(B — > A) is given by: 



C jm (B^A)= / dxdy(f> B (x,y)4>f(x,y) 
JS A 



(9) 



where Sa is the area in the plane XY enclosed by the 
hard-wall limit correspondent to A. From matrix equa- 
tions and H the generalized SM defined as: 



R A 
I B 



So S\ 
S2 S3 



I A 
R B 



can be readily obtained and is given by: 
So = -2C[B -> A)MH BA ~ 1 

Si = C{B -> A){MN + 1} 

S 2 = -2MH BA 



(10) 



(11) 



(12) 



(13) 



A. Scattering Matrix of a Single Junction. 



S, = MN 



(14) 



Mode-matching techniques can be used to obtain the 
scattering matrix (SM) for the simplest non uniform sys- 
tem: a single junction. A single junction discontinuity is 
shown in Fig. |l|, the general solution and derivatives in 
the two uniform waveguide sections are given, in matrix 
notation, by: 



<S> A {x,y,z) = I A {z) + R A (z) 
<P B (x,y,z) = T B {z) + R B (z) 



in basis A 
in basis B 



(5) 



where: 



1 = diagonal identity matrix. 



H BA = -K B 1 C T (B^ A)K A 



M = [1 — HbaC(B — > A)]' 



(15) 
(16) 
(17) 
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N = 1 + H BA C{B -> A) 



(18) 



The analysis of a single junction similar to that shown 
in Fig. [I] with the only difference that now B includes A, 
is completely analogous, equivalent equations and SMs 
are given in the next. Following two matrix equations 
are equivalent to eq. [?] and |f 



C{A -> B){I J 



R A } = I B + R B 



(19) 



P' = DI C 
R c = DR C ' 



(28) 



where D is a diagonal matrix with elements given by: 
equations 



L . With some matrix algebra from matrix 
25|and|2T 



using 



28 . one can obtain the wanted 



SM of the whole composite system: 



S A DS^DMS. 



(29) 



K A {I A -R A } = C T {A->B)K B {l B -R B } (20) 

The same definition (eq. |^) is valid for the change of basis 
matrix to change from A to B. From matrix equations 
19| and the scattering matrix, defined as stated in eq. 
10|, is given by: 



(21) 



(22) 



(23) 



(24) 



where same definitions |15[ - |18j are valid if one swaps A and 
B indexes. 



So = MN 



Si = -2MH AB 



S 2 = C(A — > B){MN + 1} 



S 3 = -2C(A -> B)MH AB - 1 



St = S A D[S B DM S A D + 1}S B 



S 2 = S B DMS A 



S 3 = S B DMS A DS B + S. 



(30) 



(31) 



(32) 



where the diagonal identity matrix has been noted as "1" 
and M is: 



B nl-l 



M = [1 - S A DS£D] 



(33) 



This same solution can be expressed in other equivalent 
ways, nevertheless, we found that this is the best formula 
to reduce computational time since it requires minimum 
number of matrix inversions: just one. 

For a greater number of junctions system, one starts 
similarly calculating each single-junction SM. Then, the 
iterative use of this combining recipe leads to the total 
SM. 



B. Scattering Matrix of many-junction structures 

The SM of a composite waveguide structure may be 
obtained by combining SMs of the corresponding junc- 
tions. The procedure to obtain the SM of a composite 
waveguide structure is described in the following. 

Let us assume two consecutive junctions joined by a 
uniform waveguide section as shown in Fig. |^. Individ- 
ual SM for both single junctions have been obtained as 
described above and are defined by: 

(£>(!?)(£) < 26 ' 

We will describe now formulae to obtain the total SM of 
the system defined by: 

(£)-(££)(£) <-> 

Since the middle connection is a uniform waveguide we 
can write: 



C. Cylindrically symmetric systems 

In the case of cylindrical symmetry, overlapping inte- 
grals get easier since can be reduced to one dimension. 
Let us detail the expressions for the simple case of plane 
hard- wall confinement. 

The symmetry assures that will be no matching of 
modes belonging to different acimutal quantum numbers 
rnHI, so we obtain a 2D equivalent problem for each m. 
The number of 2D problems that one needs to solve in- 
creases with tubes width. 

Let us consider a cylindrical uniform waveguide sec- 
tion. Electrons are confined with hard-walls at distance 
R of axis Z, the confining potential V(r) is only function 
of the distance to this axis. Thus, eq. |l]can be separated 
in cylindrical coordinates, so that the wavefunction can 
be written as: 

oo 

$ R (r,0,z)= £$£(r,0,z) (34) 

m=0 

where: 

oo 

<Z>K(r,9,z) = e* mB Y,t* l (r)F* l (z) (35) 
1=1 
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being Cm/( r ) an d Fmi( z ) to ^ e specified below. 

Due to the cylindrical symmetry, the problem can be 
considered separately for each value of m; therefore, in 
the following, we will focus in the problem for a fixed m. 

Functions (r) are the solutions of the separated ra- 
dial equation: 



1 d 



dr [ r dr^ml ( r )] 



+ {A^[E - V(r)] - [k* ml + ^jJCW = (36) 
with the boundary condition: £^?/(-R) = 

Let us consider the basis of these radial eigenfunctions 
defined between and R: B R {r) — {£mi( r )};- 12 • 
Then, the wavefunction §„\{r, 9, z) can be expressed as 
an expansion in the basis Br(t) with some coefficients 
that depend on 9 and z, written as e zme F^(z) above. In 
matrix notation this is: 



*£(r,M) 



\ 



J 



where: 



I^(e,z) = e ime I^ 



i{R)z 



\ 
) 



e R R 



(37) 



(38) 



(39) 



being k m i (R) the z component of the wavevector for the 
subband ml. 

Therefore given the (infinite) set of complex numbers 
{Rml> Imlh=i,2,- an ^ given the basis B R (r), we have 
completely described the wavefunction $^(r, 6, z). For 
any other radius R' the solution is, obviously, analogous. 

This formulation of the problem gives the correspon- 
dent SM for each m quantum number. Matrix equations 
are the same than in previous sections if we specify the 
change of basis matrix C and the transverse wavevector 
diagonal matrix K . So that for a single junction joining 
tubes of radius A and B, being B lower than A (similar 
to that represented in Fig. |l]) we have: 



CZ{B 



Jo 



iL{r)Cu(r)rdr 



[k m i(R)} 



(40) 



(41) 



being q rn i(R) given by the transversal quantized levels. 

The calculation of the relevant integral, for the simple 
case of hard-wall boundary conditions (V — Vq constant 
inside the tubes), yields: 

C%{B^A) = 

_a m h N A N B J m (am ] B/A)[J m - l { amh )~J m+l {a mh )} (42) 
2 ly mj ly mh [a mh /B]*-la m j/A] a 

In this case, radial wavefunctions are given by: 



V2 



RJm+l{ctml) 



J m {ra m i/R) 



(43) 



being a m i the zero number I of J m (x), the Bessel 
function of integer order m, and k m i(R) the z com- 
ponent of the wavevector given by: k m i(R) = 



y/ATT 2 [E-V ] 



\OLml 



/R} 2 



D. Conductance Calculations 



Within this framework, one can easily evaluate the con- 
ductance of a two terminal jSyatem on the basis of Lan- 
dauer's scattering approacnli'tj, 
formula: 



which is given by the 



G — Go T n 



(44) 



where the sum runs to all propagant modes in the input 
lead (propagant modes are defined as those with a cor- 
respondent real k, while the rest -those with imaginary 
longitudinal wavevector k- are called evanescent). Here, 
Go is the universal spinless quantum of conductance with 
the value Go = 2e 2 /h. The transmission probability for 
incidence through mode n can be easily obtained and is 
given by: 



T n = 



3 



OSa), 



2 k 



f 



(45) 



where (6*2) ,, is the (J, n) element of the S2 matrix given 
in equation pi], the sum runs to all propagant modes of 
the final lead, fcj and k n are the j and n longitudinal 
wavevectors of the final and initial tubes respectively. 



E. The truncation method 



One crucial point in the correct calculation of the SM 
of a quantum waveguide junction is the right truncation 
of two infinite series of modes representing wave functions 
at both sides of the junction. 

In the analysis of junctions, different converged re- 
sults may be obtained for different ratios of the num- 
ber of modes retained on each side. This convergence 
problem is usijiAlly referred to as relative convergence 
phenomenonllSO. Theoretical studies have found that 
the correct converged result is obtained if the ratio of 
modes on each side is taken the iS^me as the ratio of the 
corresponding waveguide widthsEH. Further calculations 
give another optimum value for this, ratio which is 1.5 
times the ratio of waveguide widthsE3. 

In our analysis we have found that one can generalize 
this optimal ratio criterium to an optimal "procedure" 
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which gives different ratios in different cases. Our ap- 
proach is based on observations of the coupling between 
different modes at waveguide junctions. 

For simplicity, let us consider the 2D axial symmetric 
hard-wall problem of a single junction like the one shown 
in the inset of Fig. |[ Assuming incidence from left to 
right and for a fixed incident mode n, the plot of the 
coupling integral C n j(w — > W) 2 versus "wide" transver- 
sal modes, j, is "peaked" at j = v such that its transverse 
wave vector satisfies q v ~ q n - Thick dots in Fig. || rep- 
resent values of this overlaping integral, the peak can be 
clearly seen at q v » q n (in this simple case q n = mr/aa 
and q v — vn/W). In a "mean field approximation"t3 
true overlaps are approximated by a uniform coupling to 
all modes within one level spacing, represented by the 
dashed line in Fig. |^. This approximation gives excellent 
total .transmission results when the difference in widths 
is big|23. We can conclude that exists a given band of 
modes of the wide side effectively coupled to the n-th 
mode of the narrow side. One needs to take into account 
all modes in this continuous interval to assure a correct 
representation of the "narrow" mode n in the "wide" ba- 
sis of transversal eigenfunctions. The truncation recipe 
that we have found is based on this observation. 

Though we have shown here a simple case, this mode 
matching scheme holds for different confining potentials 
and presumably for all of them (in particular we have 
checked some axial symmetric potentials in 2D and 3D: 
hard-wall, soft wall and harmonic. Also some non sym- 
metric 2D and 3D potentials have been checked: hard- 
wall and soft-wall) Consequently for one single junction 
similar to that in Fig. [l], the truncation procedure is 
as follows. Given transversal levels sorted in crescent 
energy: (i) Take one arbitrary number of modes N in 
the narrow side, greater or equal than the number of 
propagant ones. This implies that N can be written as 
N = N p + N e , where N p is the number of propagant 
modes and N e is the number of evanescent modes in the 
narrow tube, (ii) Take all modes of the wide side which 
are effectively connected to the ones chosen in the narrow 
side. This implies that we should use every wide-mode 
whose energy lies bellow E max given by E max = E N+1 . 
In other words, E max is the energy correspondent to the 
N + 1 level of the narrow side. 

For the truncation in the case of a greater system, in 
which the number of junctions is larger, the procedure is 
the same but choosing the narrowest of the whole system 
(i.e., the one that has the sparser transversal energy lev- 
els) as the tube to define the threshold transversal energy 
E m ax- This E max value is then used to fix the number 
of modes retained at the rest of the uniform waveguide 
sections of the system. 

Using this procedure, the number of evanescent modes 
taken in the narrowest tube N e is the only free parame- 
ter, it can be tuned for proper convergence of the results. 
In Fig. [| we can see the convergence of the conductance 
versus N e (note that the number of evanescent modes 
taken into account in the wide side could be quite larger 



than N e ). The plot has been done for the 3D "hard- 
wall" problem and is quite representative of the general 
behaviour for different confining potentials. 

The problem converges for N e equal four with an error 
of ±1%. The crucial point is that this way of making 
the truncation ensures good coupling of every mode at 
the narrow side, and the appropriate description of these 
modes is the main influence in conductance values. It 
is worth to notice here that matrices involved are rect- 
angular in general. Algebraical operations done assure 
the appropriate use of them: no numerical inversion of 
non-squared matrices has been done. One good prop- 
erty of this particular way of treating matrices involved 
is that flux is conserved with any truncation chosen; this 
implies that T n + R n = 1 even for N e = (where T n is 
the total transmission and R n is the total reflection for 
one particular propagant mode n in the initial lead). 

By the use of the matrix algebra proposed here, nu- 
merical unitarity and stability are achieved with ease. 
The SM formalism used, in which information from both 
ends of the system is taken in each step, avoids nu- 
merical roundoff errors accumulation in particular when 
combined with a suitable change of transversal basis 
as a function of the local width of the system. The 
truncation procedure explained above is the key for the 
achievement of these characteristics: it produces a self 
filtering of numerically unstable components without fur- 
ther data manipulation, hence the unitarity (to compare 
with-j'ccursion-transfer-matrix methods see for instance 



II. METALLIC NANOWIRES 

For several years now and in various experimental 
ways, the formation of ultrathin wires between two 
macroscopic pieces, of metal has been routinely achieved 
and characterized!! However, direct information about 
the atomic structure of these wires has not been available 
until recently when various experiments have achie\pd 
atomic resolution TEM images on gold nanowiresoll3. 
One of the most interesting features seen is the al- 
most ideal perfect geometry, as well as surprising heli- 
cal structuresE§E3. To the date, theoretical work has 
been focused on the study of the stability and elec- 
tronic structure for the thinnest monoatomic chainsu; 
while for the helical wires works have been focused esen- 
tially in. the prediction of the very structure for various 
metalsE.3, and-the influence of helicity in the ideal ballistic 
conductancecl In the second part of this work -by mak- 
ing use of the calculation method we have just described 
above- we approach this problem. We will point out some 
remarkable effects in the current and conductance pat- 
terns due to the ideal sharp geometries observed, these 
effects are valid in the presence or not of helicity. For a 
study on the effects of global smooth shapes ia, the con- 
ductance of atomic-scale 3D contacts see Ref.Ea. 
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A. Vortex rings in narrow-wide connections. 

A narrow-wide connection is eventually formed when 
an unsjinpprtcd gold nanowire is done by experimental 
meansQ'tj'EI This can be modelled by a simple single 
junction connection in a 3D wire similar to that sketched 
in Fig. 0. Previous theoretical work on 2D nanowires 
have pointed out the existence of a vortex structure, for 
the quantum current density, as a result of a narrow- wide 
connection. Interestingly, a change of pattern from lam- 
inar to "turbulent" 9ux is obtained for ballistic 2D free- 
electron calculationsEH. In the analogous simplest case 
for 3D wires we observe a similar behaviour and transi- 
tion, the onset of vortex structure formation is intimately 
linked to the income of the second transversal level in the 
wider part. 

In Fig. H and || we can see the charge and the current 
density for this 3D narrow-wide structure. The width of 
the narrow wire allows just a single propagating chan- 
nel while the wide part allows two channels, this is the 
minimum required to exhibit the vortex structure that 
is remarked in Fig. ^. Periodic charge and current den- 
sity structures are formed along the final wider tube, this 
can be readily explained in terms of constructive and de- 
structive interference patterns for the only two channels 
occupied in the wide lead. There are points with charge 
density zero regularly distributed at both sides of the 
charge density maximums in Fig. 0, these points are vor- 
tex centers in the quantum current map of Fig. ^. Note 
that, due to the axial symmetry, these vortex points de- 
fine rings in real space, and the flux sourrounding them 
torus surfaces. 

When the wide lead gets much wider, many transver- 
sal levels come up and the charge density structure gets 
irregular as can be seen in Fig. ^. The same explana- 
tion remains valid this time: the pattern results from 
the interference of many -in this case- longitudinal wave- 
lcnghts, the period of these patterns increases with the 
number of levels involved, eventually it gets seemingly 
non periodic as the case seen in Fig. [7] where no regular 
pattern is observed. 

These results could be important for the understand- 
ing -and therefore control- of coherent waves by tuning 
geometrical parameters. In particular the transition from 
a regular one dimensional lattice of maximums to a dis- 
ordered pattern by changing the width of the outgoing 
waveguide could be amenable of experimental verification 
(may be with electromagnetic waves rather than elec- 
trons). 



B. Resonances inside Wide-Narrow- Wide structures 

A single junction, as the one seen in previous section, 
can be regarded as a scattering center in general terms. 
Realistic wires have atomic structure and therefore many 
possible scattering centers for the conduction electrons 



passing through. The inclusion of increasing number of 
scattering centers gives rise -as we will see- to a richer 
phenomenology. In this subsection we will focus in the 
case of just two scattering centers: a wide-narrow-wide 
(WNW) junction. Two (strong) scattering centers suffice 
to predict negative differential resistance in point contact 
experiments. 

The ballistic conductance of a WNW structure, as 
the one in the,-inset of fig. ^, presents a dense reso- 
nant behaviourc3 as a function of the narrow area and 
lenght as can be seen in fig. p[ Being this a very sim- 
plistic model, it is a fair representation of what can be 
expected when we have these main ingredients: i) a cou- 
ple of strong scattering centers at a given distance, L, and 
ii) transversal confinement (characterized by the area A). 
In these cases we can expect a "wavy landscape" for the 
conductance as a function of (L, A), though may be not 
the same than the one in Fig. |[ Clearly, however, the 
variation of any of these relevant parameters can even- 
tually give any derivative. This diversity of results is 
what -we will see- is possible for the case of monoatomic 
contacts. 

Within this very simple WNW representation of 
atomic size contacts, we can easily locate the zone in 
Fig. ^ that reasonably belongs to monoatomic (or very 
few atoms) contacts for typical metals. This zone is what 
we have represented in Fig. ^| on it we have labeled three 
different points belonging to 3 different possible partic- 
ular geometries or ideal contacts. Arrows indicate the 
evolution due to an eventual stretching process of these 
few-atoms contacts. Before breaking, seems reasonable 
to think that the effective lenght—will increase and ef- 
fective area will slightly decreaseE§c£l, this is what the 
arrows indicate. Now, as can be seen, as a function of 
the starting spot, the subsequent evolution of the con- 
ductance can be very different. If we are initially located 
in point 1, the evolution will be mostly quantized (values 
very near Go) and essentially horizontal; the probability 
of this event to happen is quite high due to the great 
number of points in the plane (L, A) belonging to this 
type. However, there is a good chance of falling in a place 
like the one labeled 3, in this case -as the arrow indicates- 
we expect to observe a decreasing of G before breaking 
the contact. Finally, there is a little but non negligible 
chance for our starting point to be like the labeled 2; in 
this case, since we are located in the minimum of the 
"valley" , further stretching will make an increase in the 
conductance as indicated. 

The behaviour for case 2 is the more surprising since 
it predicts increasing of conductance while stretching 
the wire. Should this model be valid, it could be ex- 
pected to be found for given geometries and/or metals. 
In fact there is compelling experimental evidence of p3$%- 
tive differential resistance for the case of AluminumEJ'Eil, 
and theoretical explanation— in terms of atomistic den- 
sity functional calculationsE3 for the case of relatively 
wide cross-section wires. Our present free-electron result 
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provides a helpful simple view of good part of the rele- 
vant physics in these systems, a complementary piece of 
the whole picture, especially relevant for the few-atoms- 
contact case. 



reflection is naturally achieve' 



M. 



III. CONCLUSIONS 



C. Resonant Tunneling in periodic WNW structures 

Long physical wires have their thinnest expression in 
what we can call monoatomic (in width) long wires. 
The existence of -relatively long- -mono-atomic gold wires 
has beeft. exjpe jrimentally reporteda ! E3E3 and theoretically 
studiedE30cJ very recently. In a chain like this, each 
atom is a system on its own that holds localized electronic 
states in addition to those extended states responsible for 
conduction. Very simple modelling of this can be done by 
the repitition of a WNW structure. As we will see, this 
periodic repetition gives rise to a long chain of 3D "res- 
onance boxes" . With the selection of proper dimensions 
for the geometrical parameters of these boxes, localized 
states can be contained inside. The correspondent energy 
levels of this states can be an open door for resonant tun- 
neling phenomena in these 3D structures. This is proved 
here for a sufficiently long periodic wire. 

For a clear exhibition of resonant tunneling phenom- 
ena in 3D structures let us start with the standard WNW 
problem as shown in Fig. [l(]. The evaluation of G versus 
incident energy yields the expected quantized steps with 
resonance oscilations overimposed. There is a threshold 
energy at E/Ep 0.7 marking the onset of the first 
propagating level in the narrow tube, which is seen as 
the first jump in conductance. Now, if we mantain the 
same geometry but we add -in the central narrow tube- a 
wider box as shown in the inset, we can clearly appreciate 
the appearance of a resonant tunneling peak quite before 
this threshold energy. This is completely analogous to 
the obtained in ID and 2D models. 

It can be clearly prooved that perfect transparence is 
achieved at the resonant tunneling energy. Not only the 
value of transmission gets to 1 at a given E, but the pe- 
riodic repetition of this structure does not filter down 
conductivity but simply increases the number of peaks. 
There is a new peak for each 3D quantum resonant box 
we add: a splitting effect between degenerate localized 
levels is proved as can be clearly seen in Fig. 11. In this 



figure, a zoom of the resonant peak is seen for the cases 
of 1, 2 and 28 resonant boxes. 

Perfect transmission peaks are not the only conse- 
quence of interference, total reflection is also obtained as 
can be seen in Fig. [l^. A "band like" structure appears 
as a consequence of the periodicity imposed in the chain 
of 28 boxes. There are intervals of energy where trans- 
mission is nearly one and intervals where it is practically 
zero, this is in connection with conduction bands and 
energy gaps of simple one dimensional band theoriesEl 
The complete inhibition of transmission is not obtained 
in pure ID problems while for 2D and 3D problems total 



We have developed a formulation of the Scattering Ma- 
trix method to solve general three dimensional waveguide 
problems. The truncation method, based on the way the 
mode coupling is done, provides quick convergence of re- 
sults with one single parameter to be tuned: the number 
of evanescent modes in the narrowmost part of the sys- 
tem. The method has been already [Successfully used for 
the calculatioixinxiptical waveguidesEZI and electronic bal- 
listic contactsc§EjH'E£l. Examples of convergence tests 
and calculations have been also shown. In particular, we 
have studied with some detail the problem of quantum 
conductance in metallic nanowires. Various simple mod- 
els of crescent number of junctions have been studied. 
For narrow-wide connections charge density and quan- 
tum current patterns have been shown, interestingly vor- 
tex rings are obtained for this simple geometry in analogy 
to previous 2D calculations. This vortex structure starts 
forming a regular pattern and becomes apparently disor- 
dered with increasing tubes width. The calculation for 
wide-narrow-wide connections gives rise to strong reso- 
nant behaviour of the conductance as a function of width 
and length of the contact. Due to the very recent ap- 
pearance of experiments in which impressively straight 
nanowires have been imaged, seems pertinent a closer 
look at this idealized geometry model. Therefore, inter- 
pretation of WNW results, leads to the prediction of sev- 
eral behaviours of the conductance in the last plateaus of 
an experimental breaking proccess; in particular negative 
differential resistance can be expected. Further addition 
of junctions to form periodic WNW structures gives rise 
to Resonant Tunneling effects in what we can think of as 
a simple model for a long periodic nanowire. Opossite to 
ID models and analogous to 2D, not only perfect trans- 
mission is obtained for the resonant tunneling energy but 
perfect reflection appears for some finite intervals of en- 
ergy. 
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FIG. 3. Square modulus of the overlapping integral 
C„j(w — > W) versus the wide mode index j. The values 
of the widths are w/X F = 1.2, W/X F = 200.2 and the fixed 
value n taken is n — 10 (see text for more details). 




FIG. 1. Generic single junction discontinuity in a 3D quan- 
tum waveguide. Inside the hard-walls, represented here, the 
confining potential V(x, y) is arbitrary. 




FIG. 2. Building block general geometry for a composite 
waveguide structure. 
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FIG. 4. Example of convergence tests. Bottom: conduc- 
tance of a 3D cylindrically symmetric wide-narrow- wide con- 
tact versus the narrow area. Different plots can be seen for 
increasing number of evanescent modes retained in the nar- 
row tube, N e . Top: G vs N e for the worst converging point 
in the same system: A/X F fa 2.49. Parameters used in both 
plots are: L= 4 Af, A w — 4 \ 2 F . 








FIG. 5. Charge density for a NW structure in a longitu- 
dinal cut. The incident narrow tube holds a single propa- 
gant mode while two modes propagate in the wider one. The 
junction is done at z = so that, in the figure, we can see 
the charge density in the wider tube. Parameters used are 
An ■ 



0.5/Aj., Aw : 



3/A 2 F . 
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FIG. 6. Quantum current density map for the same NW 
structure in figure |^. The quantum current opens up as ar- 
rives to the wide lead, further reflection with the walls makes 
it converge to the centre. This is periodically repeated ad 
infinitum. Two sucessive magnifications give detail of one of 
the vortex ring formations. 
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FIG. 7. Charge density map 
structure than that in figure [|. 
seen analogous to that in figure I 
the axis can be seen this time. The disordered pattern is re- 
marked by contour plots at |i/>| 2 values (in axis units): 0.015 
0.02 ,0.025, 0.05 and 0.15. Parameters used are An = 0.5/A| 
Aw = 200/A|. 




FIG. 8. Conductance for a 3D WNW cylindrically sym- 
metric contact (after Ref. [28]). 
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FIG. 10. Conductance vs incident Energy for WNW (dark 
line) and a Resonant Tunneling single box structure. Both 
geometries are shown in the inset. A clear total transmission 
peak appears before the onset of the first propagating level 
inthe narrow tube. As well, the resonant box induces total 
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FIG. 11. Resonant Tunneling peak magnified for 1,2 and 
28 boxes (see text). 
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